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Magnetic thin films exliibit a strong variation in properties depending on their degree of disorder. 
Recent coherent x-ray speckle experiments on magnetic films have measured the loss of correlation 
between configurations at opposite fields and at the same field, upon repeated field cycling. We 
perform finite temperature numerical simulations on these systems that provide a comprehensive 
explanation for the experimental results. The simulations demonstrate, in accordance with exper- 
iments, that the memory of configurations increases with film disorder. We find that non-trivial 
microscopic differences exist between the zero field spin configuration obtained by starting from a 
large positive field and the zero field configuration starting at a large negative field. This seem- 
ingly paradoxical behavior is due to the nature of the vector spin dynamics and is also seen in the 
experiments. For low disorder, there is an instability which causes the spontaneous growth of line- 
like domains at a critical field, also in accord with experiments. It is this unstable growth, which is 
highly sensitive to thermal noise, that is responsible for the small correlation between patterns under 
repeated cycling. The domain patterns, hysteresis loops, and memory properties of our simulated 
systems match remarkably well with the real experimental systems. 

PACS numbers: 



I. INTRODUCTION 

Many magnetic systems have a memory of their past 
configmations. This history, manifest in the hystere- 
sis loop of a magnet, has many fascinating features, 
for example Barkhausen noise which demonstrates that 
this history has intricate behavior at a small scale, with 
avalanches occurring in the same order in a highly re- 
producible manner [lllllllllllEllli. From a the- 
oretical perspective, a certain class of Ising models has 
been proved by Sethna et al. Q to exhibit perfect return 
point memory 27] And of course, hysteresis is the prin- 
ciple that makes magnetic storage devices possible [lOl |. 

Primarily due to the technological importance of mag- 
netic systems, a vast number of experiments have been 
done to observe magnetic memory and hysteresis phe- 
nomena ^ comparably large number of theories 
have been made to explain these experiments with vary- 
ing de gree s of success. A recent experiment by Pierce 
et al. pdl Il2| measured the memory properties of mag- 
netic multilayer thin films. They observed an effect that 
at first sight appears paradoxical, involving fundamen- 
tal issues of symmetry in these systems. In this paper, 
through numerical simulations, we attempt to provide an 
explanation for this effect. 

The experiment by Pierce et al. used the powerful new 
tool of x-ray speckle metrology to measure the covariance 
of a domain configuration at one field with the configura- 
tion at another field. This measured value is one way to 
quantify the amount of "memory" possessed by the mate- 
rial. Measurements were done for Co/Pt multilayer sam- 
ples with varying amounts of disorder; the study found 
that samples with greater disorder had higher "mem- 
ory" than the ordered samples which showed no signif- 
icant amount of memory. The domain patterns ranged 




FIG. 1: Magnetization vs. external field for a cartoon hys- 
teresis system. Return point memory is defined as the covari- 
ance between any point with the same point after an integer 
number of complete cycles. Complementary points are points 
lying a half-integer number of field cycles away. Each pair of 
points (filled, shaded, and unfilled) is a pair of complementary 
points. 



from labyrinthine mazes for low disorder samples to ones 
without any noticeable structure for the highest disorder 
samples. Hysteresis loops for different samples had dra- 
matically different features depending on the amount of 
disorder, with steep cliffs (large changes in the magneti- 
zation at constant field) existing in the hysteresis loops 
of the low disorder samples. 

In addition to the memory dependence on disorder, 
an unexpected finding was the difference between what 
Pierce et al. called "return point memory" (RPM) and 
"complementary point memory" (CPM). RPM was de- 
fined by them as the covariance, see below, between the 
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configuration at a certain field point with the configura- 
tion at the same field point after an integer number of 
complete cycles around the major hysteresis loop. CPM 
is defined as the covariance between the configuration at 
a certain field point with the configuration at a half inte- 
ger number of field cycles away. Figure ^ shows examples 
of "return" and "complementary" points. The covariance 
that we will use in order to define RPM and CPM is 

cov{a, b) = (sa(r) • Sh(r)) - (sa(r)) • (sb(r)), (1) 

where a and b refer to legs on the hysteresis loop and 
the average is over all space (spins). The normalized 
covariance is defined by 

p=_^^±l^^. (2) 

■\/ cow(a, a)cov{b, b) 

At a certain external field B^, RPM is defined as 
p{Be, a; Be,b), where a and b are legs going in the same 
direction (both ascending or both descending) , and CPM 
is defined as —p{Be, a; —Be, b), where a and b are legs go- 
ing in opposite directions. Perfect memory would lead to 
both RPM and CPM equal to unity. A covariance be- 
tween two independent systems would approach zero for 
both RPM and CPM. 

The speckle experiment measured intensities from x- 
ray scattering. Because scattering experiments probe 
samples in Fourier space, the actual measurements did 
not measure the covariance as defined by Eq.lQ) but in- 
stead by 

covu{a,b) = (|.J(k)n.^(k)r)-(|<(k)r)(|.,^(k)p), (3) 

where Sa(k) is the z component of the Fourier transform 
of the spins on leg a and the the average is over wave 
vectors |23|. Using this Fourier space covariance in the 
definition of RPM and CPM, Pierce et al. measured the 
amount of memory for samples of Co/Pt multilayer thin 
films with a wide range of disorder. RPM and CPM 
are both negligible for samples with low disorder and are 
quite significant (> 0.6 at the coercive field) for highly 
disordered samples. The onset of memory occurs quite 
abruptly as disorder is increased. Furthermore, for sam- 
ples with significant memory, RPM is noticeably larger 
than CPM for all fields measured. In this paper we will 
use the real space covariance in the definition of RPM 
and CPM. 

A viable model of the Co/Pt multilayer films must at 
the least be able to explain two experimental results: 

1. RPM, CPM ~ for systems with low disorder 

2. 1 > RPM > CPM for systems with high disorder 

In this paper, we provide a possible explanation for these 
experimental results through finite temperature numeri- 
cal simulations of classical vector spins. The finite tem- 
perature destroys perfect memory and affects both RPM 
and CPM. We will see that the nature of the domains 



determine how much the temperature affects the covari- 
ance; highly ordered systems are more susceptible to tem- 
perature effects than highly disordered systems, and do 
explain result ^ But thermal noise does not discriminate 
between the ascending leg and the descending leg. A 
satisfactory theory must also show that RPM is greater 
than CPM. 

Condition [2 is difficult to understand intuitively. Con- 
sider a system at low temperature that starts out fully 
saturated from being in a large positive external field 
perpendicular to the film, and then that field is brought 
down adiabatically to zero. At remanence the system 
will be in a state with domains pointing in different di- 
rections. Now repeat the same procedure but starting 
with a fully saturating field pointing in the opposite di- 
rection. At remanence we now expect the final configura- 
tion to be the same apart from a change in sign s — s. 
The fact that the CPM value is less than the RPM value 
contradicts this. The configurations although somewhat 
correlated, are different. 

One possible explanation of this can be devised by in- 
troducing a term in the Hamiltonian that is not invari- 
ant under spin and external field reversal s ^ — s and 
Bf, — > —Bf,. Terms in this class automatically introduce a 
difference between the ascending and descending legs. A 
system with such a Hamiltonian would exhibit a hystere- 
sis loop that is not symmetric under the same operations 
even at zero temperature. Using a scalar 0* model with 
a random field term in the Hamiltonian, Jagla [13 . 
was able to satisfy both conditions. But the physical ori- 
gin, and existence, of these random fields is unclear. Re- 
cently, with some modifications to his 0'* model, includ- 
ing replacing the random fields with random anisotropy, 
Jagla [iJI was able to produce domain patterns and hys- 
teresis loops that resemble the experimental patterns and 
loops remarkably well. However without a random field 
term in the Hamiltonian, the memory conditions could 
not be satisfied. 

We provide a possibly more fundamental mechanism 
that satisfies condition |2 The vector dynamics breaks 
the spin and field reversal symmetry, thereby reducing 
CPM and not RPM. This mechanism does not require 
any new terms in the Hamiltonian of the class men- 
tioned in the previous paragraph. In our earlier work |l5l | , 
we showed how the vector dynamics, governed by the 
Landau-Lifshitz-Gilbert (LLC) equation can give 

rise to non-complementary hysteresis loops for a system 
of nanomagnetic pillars. The crucial role of vector dy- 
namics for explaining experimental results underlines the 
inadequacy of scalar models, even when the system is 
highly anisotropic. We have not excluded the possibility 
that the real explanation is a combination of the vector 
dynamics and the random fields. Furthermore it is pos- 
sible that the experiments result from some other effects 
that have not been considered so far. In this paper, we 
use numerical simulations to demonstrate the plausibility 
of this vector mechanism. 

In the next section of this paper, the LLC equation is 
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introduced and the various terms in the Hamiltonian and 
their origin will be described. Details of the numerics are 
provided in Sec. Ill and the types of domain patterns and 
hysteresis loops from these simulations are presented in 
Sec. IV. Sec. V contains the covariance results. 



II. CLASSICAL SPIN DYNAMICS AND THE 
MODEL HAMILTONIAN 

The Landau-Lifshitz-Gilbert equation of motion 
is the simplest dynamic equation for classical spins that 
contains a reactive term and a dissipative term: 



ds ,^ ^, 

— = -s X (B -7s X B), 



(4) 



where s is a microscopic classical spin, B is the local ef- 
fective field, and 7 is a damping coefficient. The effective 
field is B = —dTi/ds + where 7i is the Hamiltonian 
and C represents the effect of thermal noise. The overall 
constant in front of the right hand side of Eq.Q is set 
to unity. 

Both terms on the right hand side of the LLG equa- 
tion are necessary to adequately describe the motion of 
classical spins. The double cross product in the damping 
term is necessary in order to maintain a constant mag- 
nitude for the spins. The relaxation time is inversely 
related to the damping coefhcient. The reactive term de- 
scribes precession of a spin about its local effective field 
and has a quantum mechanical origin; the commutation 
relations of angular momentum variables give rise to the 
single cross product of the reactive term. The commu- 
tation relation between spin variables is odd under spin 
reversal. Later, we will show how precession is crucial to 
adequately describe the experimental system using our 
Hamiltonian. Because scalar models cannot have preces- 
sion, the mechanism that we suggest is not possible in 
scalar theories. 

The LLG equation can be applied to microscopic spins 
as well as coarse grained spins. Every spin variable in 
the numerics represents a block spin of the multilayer film 
and the evolution of all spins is calculated by numerically 
integrating the LLG equation. Henceforth in this paper, 
a "spin" corresponds to a block spin variable. The coef- 
ficients of the Hamiltonian, the damping coefficient, and 
thermal noise are also those associated with the coarse 
grained variables. The Hamiltonian has four terms: self 
anisotropy energy, local ferromagnetic interactions, long 
range dipole-dipole interactions, and the energy for the 
interactions with the external field. 

The Co/Pt multilayer film is a perpendicular 
anisotropic material. We define this perpendicular axis 
as the z axis. The origin of the perpendicular easy axis is 
the layered structure of the material. Any real material 
will have imperfections in the layering, and this imper- 
fection in the layering is the "disorder" mentioned many 
times in the previous section. Samples with low disor- 
der show very well defined planes separating the layers 



of the two elements, whereas samples with high disorder 
have very rough interfaces between the Co and Pt. Due 
to this disorder in the layering, each spin has a different 
easy axis, n,;. The anisotropy energy term is described 

by 



'Ham = -a^{si ■ hi)'^ 



(5) 



where a is a model parameter. The anisotropy must be 
an even function due to the symmetry in ±11^. Aside 
from being even, the functional form of Eq.(|SJl can be 
quite complicated in general, but if a power expansion is 
possible, Eq.Q would be the leading order term. 

The next term in the Hamiltonian describes the local 
ferromagnetic coupling between neighboring spins. 



(6) 
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where J is the ferromagnetic coupling constant. We con- 
sider only the J > case. Even though we attempt to 
model a continuum system, a grid is necessary for the 
numerics. In order to minimize artificial effects due to 
the grid, this real space ferromagnetic energy term is re- 
placed by a term in Fourier space of the form 



He 
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(7) 



Eq.((7| is the leading order term of the elastic energy in 
Fourier space which aligns spins locally similar to Tij. 
Disorder can be introduced in the elastic constant J as 
well. For small disorder of this type, the results presented 
in this paper do not change qualitatively. We eliminate 
this extra disorder parameter and have a uniform elastic 
constant for all spins. 

In addition to the local interactions, the spins also in- 
teract via a long range dipole-dipole energy. The form of 
this energy is the usual classical expression 



dip 



-w J2 



(8) 



iwhere is the displacement vector between spins i and 
J, eij is the unit vector along this direction, and w is 
another model parameter. When the spins are point- 
ing predominantly in the ±z directions, the dipolar fields 
tend to anti-align the spins. This competition with the 
ferromagnetic interaction produces many of the domain 
features. Interestingly, when the spins are in-plane, the 
dipole-dipole and local interactions can become cooper- 
ative (depending on r^j). 

The form of the dipole energy expressed by Eq.® is 
correct for point dipoles, but must be corrected for the 
block spins of interest. The small but finite thickness of 
the layers does not significantly change the long range 
behavior of Eq. (|HJ| , but it does eliminate the divergence 
at small separations . This cutoff is implemented in 
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Fourier space which is a more computationally efficient 
basis to calculate long range dipolar fields. Multiplying a 
gaussian of the form exp(— to the Fourier trans- 
form of Eq. (|SJ) effectively removes the divergence at large 
wave vectors while retaining the original form for small 
wave vectors. The model parameter, d, is the length scale 
of the cutoff where the finite thickness of the sample af- 
fects the dipolar interaction. 

The last term in the Hamiltonian is the energy from 
the interaction of the spins with the external field which 
we take to be uniform in the z direction, 

ne.t = -BeY,sl. (9) 

i 

Major hysteresis loops are simulated by cycling the ex- 
ternal field from a large positive value (past saturation) 
to a large negative value and back again. To ensure to- 
tal saturation we reset the value of the spin variables to 
point completely along the field direction at the start of 
each leg in the hysteresis loop. 

In summary, the full model Hamiltonian has four 
terms: anisotropy, ferromagnetic coupling, dipole-dipole 
interaction, and the external field term as described by 
Eqs.HaH), 

= Hani +^£1 +'Hdip + 'Hext- (10) 

The model Hamiltonian is bilinear in the spins, and ex- 
ternal field, Be, therefore it is symmetric under — 
and Bf. — » —B^. However, as mentioned above, this sym- 
metry does not exist for the equation of motion. Un- 
der these two operations, the left hand side of the LLG 
equation, Eq.Q), and the dissipative term reverse signs, 
whereas the reactive term stays unchanged. With the 
symmetric model Hamiltonian and using purely relax- 
ational dynamics, in some sense there would be no way 
to differentiate from going down the descending leg or up 
the ascending leg of the hysteresis loop. The difference 
in sign change between the terms of the LLG equation 
under spin and field reversal is the mechanism by which 
RPM > CPM. When precession is turned off, the covari- 
ance results presented in Sec. V are no longer valid and 
we verified that as expected, RPM CPM when the tem- 
perature is non-zero and RPM — CPM = 1 identically 
when there is no thermal noise. 

This explanation for why RPM > CPM may become 
invalid when the relaxation time is much smaller than 
the precessional period. There is no reason to believe 
that this is the case for multilayer thin films. In fact, 
7, the damping coefficient which is a measure of the rel- 
ative importance between damping and precession has 
been measured for NiFe thin films and found to be ap- 
proximately 0.01 This indicates that the relaxation 
time is much larger than the precessional period, thus 
precession is substantial for some systems. The damping 
is found to be enhanced to ~ 1 for CoCr/Pt multilayer 
films 'IQ]. For the numerical simulations, we have set 7 
to unity. 



III. NUMERICS 

We simulate the multilayer thin film by a two dimen- 
sional lattice of block vector spins of unit length. Most 
of the model parameters are noted in the previous sec- 
tion: the relative strengths of the coupling constants in 
the Hamiltonian, a, J, and w, the dipolar cutoff length 
d, and the damping coefficient 7. There are two more 
model parameters, the temperature, T, and A, a parame- 
ter that controls the variation in the easy axes, rii. Each 
fii is assigned randomly for every spin. The stacking of 
the layers in the z direction is the physical origin of the 
anisotropy, therefore the n^'s are weighted in this direc- 
tion. This weighting factor. A, is inversely related to the 
amount of disorder. 

The large number of parameters may at first seem 
hopeless from the point of view of prediction because in 
a seven dimensional parameter space almost any curve 
can likely be fit. This clearly diminishes the predictive 
power of a model. All is not lost, because the behav- 
ior described in this paper is quite robust in many of 
these parameters. Furthermore there is agreement, at 
least at the qualitative level, for all the quantities mea- 
sured at a fixed value of parameters: the shape of the 
hysteresis loops as a function of disorder, the evolution 
at low disorder of patterns in the films which have also 
been seen experimentally, and the RPM and CPM be- 
havior as a function of disorder. As we mentioned at 
the beginning of the last section, 7 has been measured 
in a certain multilayer system to be ~ 1 and we 

have set it equal to unity for all simulations. A smaller 
7 would almost certainly enhance our results. The dipo- 
lar cutoff, d, determines a length scale of the domains, 
and all other parameters are considered with respect to 
this length scale. Without loss of generality, we have set 
d = A lattice spacings. Below, we will report results for 
different temperature and disorder, therefore T and A are 
not fixed model parameters. That leaves us with a three 
dimensional parameter space and the task does not seem 
so daunting as before. 

One may ask that since we are simulating a known ex- 
perimental system, why are we not using measured values 
for these quantities in our numerics? First, these are not 
known for these disordered Co/Pt films, and even if this 
could be done, the simulations are of coarse grained spins 
and not microscopic spins therefore all coefficients, for 
example, temperature, magnetic field, disorder strength, 
would change in a highly nontrivial manner. Lattices of 
sizes 128 x 128 and 256 x 256 spins are used in the simula- 
tions, therefore each block spin represents roughly forty 
thousand real microscopic spins. Nevertheless, even with 
only 128^ spins, we are able to observe behavior quite 
similar to what is seen in the experiments. 

Initially, the easy axes, fii are independently and ran- 
domly assigned for each spin. The z component is 
weighted by A which sets the amount of disorder; disorder 
is small when A is large and is increased by decreasing A. 
The system is initially saturated in the positive z direc- 



tion. Starting from a large positive external field, 
the field is adiabatically decreased to a large negative 
field, —5™"^ and the spins are again saturated. From 
there, the field is adiabatically increased back to i?^"*^ 
to finish one complete field cycle. 

At zero temperature, adiabatic field cycling is straight- 
forward: after every field step, the spins evolve in time 
until the configuration reaches a local minimum in the 
energy, then the field is changed by another step. With 
thermal noise, the adiabatic condition cannot be as strin- 
gent as it is at absolute zero. The adiabatic condition 
becomes nontrivial because thermal noise is constantly 
buffeting the spins. The systems of interest have a large 
anisotropy with easy axes weighted in the z direction, 
therefore the energy minima of the spins are predomi- 
nantly in the ±z-direction. If the temperature is not too 
large, it would be unlikely for thermal noise to be suffi- 
cient to knock the spin over the energy barrier to induce 
a spin-flip. 

We use the fact that thermally activated spin-flips are 
uncommon to construct an adiabatic condition. The 
spins are allowed to vary their orientation from the ver- 
tical, but if the z component of any spin changes from 
greater than a large threshold value to less than the neg- 
ative of that same value (or vice versa), then the system 
has not yet satisfied the adiabatic condition. If spin-flip 
events do occur, the spins must be evolved further with- 
out changing the field. This threshold value is set to 
0.75 for the results in this paper. Qualitatively similar 
results are obtained for slightly different thresholds. In 
addition to not allowing large spin flips for the adiabatic 
condition, the change of the total magnetization in the z 
direction, AM, must also not change substantially. The 
condition used is AAf < 0.1. 

The time evolution of the spins follows the LLG equa- 
tion. This is done by calculating the effective field for 
each spin and then integrating the equations using the 
fourth order Runge-Kutta algorithm. Due to the long 
range nature of the dipolar interaction, all the spins are 
coupled together and a real space numerical integration 
is O(iV^). Using a FFT algorithm, this is reduced to 
0(N\a.l>S\ Periodic boundary conditions are used to 
accommodate the Fourier method, but the results here 
should not differ significantly with different boundary 
conditions. 

The configurations of spins are stored at certain fields 
and RPM and CPM as defined in the first section are cal- 
culated between different legs. For efficiency, many legs 
are run in parallel. In the next section, the domain pat- 
terns and hysteresis loops resulting from these numerics 
are discussed. 



IV. DOMAIN PATTERNS AND HYSTERESIS 
LOOPS 

Using the numerics described in the previous section, 
qualitatively different configurations and hysteresis loops 
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FIG. 2: Domain growth for a 256^ system with low disor- 
der, A = 1000 and w = 0.15. The z component of the 
spins are shown with light shades representing "up" and dark 
shades representing "down." The first five panels are at con- 
stant field and show snake-like domains growing in time. The 
labyrinthine patterns resemble the domain patterns of the low 
disorder samples from the speckle experiments. Eventually, 
the length of the domains becomes comparable to the system 
size and require a lowering of the external field to increase in 
size. The last panel is at a lower external field than the other 



arise depending on the parameters described in Sec. II. 

The important parameters are the relative strengths be- 
tween the interaction terms in the Hamiltonian (a, J, 
and w), and the amount of disorder (characterized by A). 
The coarse graining introduces difficulties in calculating 
these parameters from the microscopic interactions. Fur- 
thermore, even a first principles calculation of the micro- 
scopic interactions is highly nontrivial. For example, the 
anisotropy of the spins is believed to be due to a quan- 
tum exchange mechanism between layers that is not eas- 
ily calculated or understood. Nevertheless, we can pro- 
ceed by empirically searching through parameter space 
to find parameters in which the experimental results are 
observed. 
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In the region of parameter space in which the domain 
patterns and hysteresis loops have properties seen in the 
real experimental samples, the ferromagnetic coupling J 
and the anisotropy strength a are approximately equal, 
whereas the dipole-dipole strength w is roughly an order 
of magnitude weaker. The competitive nature of these 
interactions give rise to interesting domain patterns and 
hysteresis loops 0, A strong a tries to keep all 

spins pointing out of plane due to the weighting in the z 
direction. The positive J tries to maintain a local align- 
ment of spins, thus large domains are favored by this 
energy term. Though weaker in magnitude per spin-spin 
interaction, the dipolar force is long range with an accu- 
mulated effect that competes significantly with the local 
couplings. 

The amount of disorder in a system is determined by 
A, with a large A signifying low disorder. When A = 1, 
the anisotropy energy is spherically symmetric (after av- 
eraging over many spins) and the system is no longer a 
perpendicular material. In addition to lowering A, dis- 
order also weakens the dipolar strength. Disorder in the 
real experimental systems exists in the interface between 
the Co and Pt layers. When disorder is large, the in- 
terface is highly non-planar, thus effectively suppressing 
the dipolar interaction after coarse graining. Though this 
disorder in the interface would most likely alter the rela- 
tive strengths of J and a, we remove this degree of free- 
dom and set J = 0.85 and a = 0.875 for all systems. 

The low disorder systems have the most interesting 
domain patterns with labyrinthine mazes at remanence. 
The systems with the lowest amount of disorder studied 
have A = 1000 and w = 0.15. For this value of A, the 
easy axes is almost parallel to z. Starting at positive 
saturation, the large J and a keep the spins aligned es- 
sentially vertically. But the dipole-dipole interaction is 
highly dissatisfied in this configuration. Because of this 
dissatisfaction, at a relatively high positive field, one or 
more small patches will flip. 

Once a small "down" domain has formed in a sea of up 
spins, the local ferromagnetic interaction causes instabil- 
ities to occur at the boundary of the domain. Due to the 
ordered environment around the domain, the instabilities 
grow into lines. These configurations are long serpentine 
configurations in accordance with the recent experiments 
on Co/Pt films and also earlier experimental work [20l |. 
The growth of the snakes is spontaneous, occurring at 
constant external field. Figure |21 shows the formation 
and growth of the snakes. In regions of high curvature, 
occasionally side branches grow and look very similar to 
experimental patterns. Due to the dipolar force, the do- 
mains behave as if they repel each other. This repulsion 
causes the growth to halt once the length of the snakes 
is comparable to the system size and the total area of 
negative spins fills a finite fraction of the sample. At this 
point, in order for more spins to flip, the field must be 
lowered. 

The initial circular instability and subsequent growth 
occurs for a wide range of temperatures and therefore ap- 



pears to be predominantly different from what one would 
expect if this was nucleation. If this was nucleation, near 
the critical field, domains would appear and disappear 
with a temperature dependent probability. Our simula- 
tions show that when a domain flips initially, it remains 
in that state above the critical field and never disappears. 
At the critical field a domain will quickly become unsta- 
ble, even at zero temperature, and grows into serpentine 
patterns even with no change in applied field. The ini- 
tial instability often occurs at the same location upon 
repeated field cycling which further indicates that this is 
not nucleation. 

The growth of these snake-like domains is similar to 
what is seen in ferrofluids [TtL l2ll . In that case how- 
ever volume of fluid is conserved which means that the 
volume of fluid cannot spontaneously increase with ex- 
ternal field. However the general shape of domains is 
similar, particularly their long serpentine shape ending 
with a wider head as shown by figure [3 

The growth at constant field introduces a cliff in the 
hysteresis loop; the growth of the snakes reduces the mag- 
netization abruptly. Hysteresis curves of this type are 
shown in the first panel of figure 21 The end of this 
large avalanche signifies the point at which the snakes 
have taken up the available growth environment and are 
intertwined. After this point, the length of the snakes 
can no longer increase significantly but the width of the 
snakes does gradually increase by lowering the external 
field. Eventually, the field overcomes the repulsion and 
the snakes link and form a labyrinthine maze. Finally, 
the field will have a large negative value that saturates 
the spin in the negative z direction. 

The systems with the highest amount of disorder stud- 
ied have A = 2 and w — 0.05. With such a low weight- 
ing, the easy axes have large components in the plane 
that vary greatly from one spin to the next. Because the 
hi 's no longer have a large probability of pointing in the z 
direction and the anisotropy energy is strong, the satura- 
tion magnetization is reduced. In other words, the spins 
tend to point along their easy axes which have larger in- 
plane components than the low disorder systems. This 
reduction of the saturation magnetization as disorder is 
increased is also seen in the experiments. Though the 
spins are more in-plane, the external field still starts 
at a large positive value, therefore the z component of 
the spins start positive. The weaker dipole coefficient 
requires more reduction of the field before the dipolar 
dissatisfaction causes a patch to flip, thus the remanent 
magnetization increases for more disordered systems. 

Disorder causes the domains to have a rougher bound- 
ary. This rough edge is due to the more spherical sym- 
metric easy axes of the spins at and near the edge. Fur- 
thermore, the highly irregular easy axes landscape pre- 
vents long snake-like domains from forming. Therefore, 
domains cannot grow significantly unless the external 
field is lowered. The stunting of the domain growth can 
be seen by a lack of any large avalanches seen in the 
hysteresis loop. 
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FIG. 3: Spin configurations at the coercive field for systems 
witli different amounts of disorder. Tlie systems witli tlie 
least amount of disorder are in the upper left panel and 
greatest amount are in the lower right. The disorder pa- 
rameters A are 1000, 11, 4.1, 3 and the dipole strengths w are 
0.15,0.105,0.08,0.06. 




FIG. 4; Hysteresis loops for systems with different amounts 
of disorder; the vertical axes are the magnetizations and the 
horizontal axes are the external fields. The panels correspond 
to the same system as in figure The cliff in the hysteresis 
curve vanishes as disorder is increased. 



Figure 01 shows configurations near the coercive field 
for systems with increasing disorder. Figure 0] shows the 
corresponding hysteresis loops for these systems. The 
configurations and hysteresis pfots are for systems with 
A = 1000, 11, 4.1, 3 and w = 0.15, 0.105, 0.08, 0.06 start- 
ing from the upper left to the lower right panels. This se- 




External Field 



FIG. 5: Hysteresis loops at different temperatures for a low 
disorder system (top, A = 1000, = 0.15) and a high disor- 
der system (bottom, A = 4.1, ui = 0.08). For both systems, 
thermal noise narrows the hystersis loop. For the low disorder 
system, on the descending leg thermal noise causes the cliff 
to occur at a higher field. 



quence of hysteresis plots shows how the height and slope 
of the cliff decreases as disorder increases. At an inter- 
mediate amount of disorder, snakes do grow but not to 
the extent of the entire system size. The domain growth 
is pinned by the disorder and growth cannot occur un- 
til the field is lowered. Eventually, beyond A = 4.1 and 
w = 0.08, there is no noticeable cliff in the hysteresis 
curve. A system with a lower w remains saturated for 
a larger range of external field, therefore the remanent 
magnetization and the width of the hysteresis loop in- 
creases as w decreases. 

The general features of the domain patterns and hys- 
teresis loops described above are dependent upon the 
amount of disorder but are independent of temperature 
for a wide range of temperature. However, thermal noise 
does affect the field in which the domains initially flip. 
This effectively narrows the hysteresis loops as shown in 
figure |31 Another study of Co-Pt multilayers shows this 
effect as well 0|. Though thermal noise does have an 
effect on the hysteresis loop, domain creation and growth 



is not predominantly nucleation. 

One can understand this by examining the low disorder 
system. Starting from positive saturation, at a critical 
field, circular initial domains and serpentine growth oc- 
curs for all temperatures. At zero temperature nucleation 
cannot occur, therefore the zero temperature critical field 
is analogous to the spinodal point. When thermal noise 
is present, it can cause domains to flip before the spin- 
odal field is reached, thereby nucleating domains which 
then grow. Hence at finite temperature, the cliff occurs 
at a field larger than the zero temperature critical field 
(when all spins are initial saturated in the positive direc- 
tion). But from figure [S] we see that this is not a strong 
effect. Nucleation only slightly changes the field in which 
the cliff occurs and is not the dominant phenomenon for 
domain formation and growth. 

The qualitative features described in this section are 
observed in the experimental samples. In summary, low 
disorder systems have labyrinthine domain patterns and 
hysteresis loops with steep cliffs and high disorder sys- 
tems have irregular patchy domains and more "standard" 
hysteresis loops. 



V. COVARIANCE RESULTS 

As seen in the previous section, the amount of dis- 
order and dipole strength dictates the type of domain 
patterns and hysteresis loops. In this section, we discuss 
how memory properties, specifically RPM and CPM, are 
determined by the amount of disorder. Though thermal 
noise did not have a large effect on the hysteresis loops 
and domain pattern, memory properties are highly sen- 
sitive to temperature. 

At zero temperature, the snakes seen in the domain 
patterns of the low disorder systems meander due to the 
small differences in one direction versus another. These 
slight differences are mainly due to the randomness in the 
easy axes of the spins. Without thermal noise, every as- 
cending leg will be identical to the next; a complete field 
cycle that saturates the magnet returns the system to ex- 
actly the same state. In other words, RPM is identically 
equal to unity. 

Because the meandering is due to very slight differ- 
ences in the local environment of the snake, thermal noise 
could alter the domain patterns dramatically. Even when 
the temperature is low, the domain configurations can be 
drastically different between points separated by an inte- 
gral number of complete field cycles. In the second row 
of figure ini two configurations of a low disorder system 
are so different that one could not tell that both panels 
are the same system separated by one field cycle. These 
two configurations both appear, at least visually, to have 
domains that are uncorrelated both in their positions and 
their shapes. 

Of course high disorder systems are also affected by 
thermal noise. But because these systems are much more 
heterogeneous, a small amount of noise does not cause 




FIG. 6: Comparisons of the spin configurations after one com- 
plete cycle at the coercive field. The configurations on the 
top row are from a system with high disorder (A = 4.1 and 
w — 0.08). Notice the domains have essentially the same 
shape and are in the same positions. The bottom two panels 
are configurations from a system with low disorder (A — 1000 
and w = 0.15). There does not appear to be any obvious 
correlation between the domains for the low disorder system. 
The temperature is 10~^. 



such drastic differences. The lack of domain growth at 
constant field is due to the orientational disorder. This 
disorder also ensures that even for different realizations of 
thermal noise, the domains form in essentially the same 
locations with the same sizes. In other words, the rel- 
atively more reproducible growth of domains is another 
consequence of pinning by disorder. Therefore the differ- 
ences in configurations after a complete field cycle, shown 
in the top row of figure IHl is not as drastic as was seen in 
the low disorder case. 

Whereas thermal noise adequately explains return 
point differences, explanations of the differences between 
complementary points also require a discussion of the dy- 
namics. When the precession term of the LLG equation 
is removed, the behavior of the complementary branch 
is identical to the return branch (with s — > — s and 
Be — *■ —Be). But when precession is present, even at 
zero temperature, a configuration at field B along the as- 
cending branch is not identical (after spin reversal) to the 
configuration at field —B along the descending branch. 

As mentioned in the introduction, RPM and CPM 
quantify the amount of correlation between configura- 
tions. Because each system is saturated after every leg, 
RPM and CPM are independent of the number of inter- 
mediate field cycles. At very low temperatures, the co- 
variance values for different legs essentially overlap. Fig- 



FIG. 7: RPM and CPM vs. external field for a A = 4.1 
and w — 0.08 system at a very low temperature of 10~^. 
The covariance values for four "return" legs (dashed lines) 
and four "complementary" legs (solid lines) are plotted. All 
"return" legs have essentially the same covariance. Similarly 
all "complementary" legs have the same covariance. Because 
the system reaches saturation at the end of every leg, RPM 
and CPM are independent of the number of intermediate field 
cycles. 



ure[7|shows RPM and CPM plots for a highly disordered 
system with A = 4.1 and w — 0.08 at a very low tempera- 
ture T = 10"^. RPM is greater than CPM for most field 
points and the different RPM legs lie very close together 
as do the different CPM legs. This feature is seen in the 
experiments. 

The effects of temperature on RPM and CPM are 
shown for a system with A = 4.1 and w = 0.08 in fig- 
ure |S| There are much larger fluctuations compared to 
the system in figure|7|due to the higher temperature. As 
expected, these fluctuations are reduced as the system 
size is increased. The average is done over many legs and 
the root mean square deviation is shown. These figures 
clearly show that at low temperature RPM is close to 
unity whereas CPM is significantly lower. As tempera- 
ture increases, RPM decreases as discussed. Curiously, 
for low temperatures, within errors, CPM stays constant. 
Up to a temperature of ~ 10~^, the dynamics appear to 
be dominant over temperature in terms of complemen- 
tary point memory. This result emphasizes the impor- 
tance of the vector dynamics. Furthermore, in the region 
with RPM, CPM < 1 the amount of memory decreases 
with field. This decrease in memory with field is seen in 
the experiments as well. 

For a temperature of lO"**, figure El shows the aver- 
age RPM and CPM for systems with different amounts 
of disorder. The low disorder system has no significant 
amount of memory near the coercive field. Furthermore, 
a spike exists at the critical field where the instabilities 
initially form. This spike reveals the fact that the ini- 
tial domain(s) forms at the same position after repeated 
field cycling, hence the large covariance. RPM and CPM 
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FIG. 8: RPM and CPM vs. external field for a highly dis- 
ordered system (A = 4.1 and w = 0.08) at different tem- 
peratures. At very large negative fields, temperature is in- 
versely related to the height of the lines. The temperatures 
are T = 0, 10"^ 10"", 10"^. RPM and CPM values shown are 
the mean values for multiple field cycles with RMS error bars. 
The top panel shows how RPM is exactly equal to unity for 
zero temperature and decreases as temperature increases as 
expected. Unexpectedly, as shown in the bottom panel, CPM 
does not always decrease as temperature increases. In fact, 
except for the T — 10~^ case, CPM is essentially constant for 
all other temperatures within errors. Because CPM is similar 
between the zero and small temperature cases, thermal noise 
is not the dominant cause of loss of complementary memory. 



drop off sharply away from this critical field due to the 
thermal sensitivity of the domains. For A < 4.1 and 
w < 0.08 both RPM and CPM are noticeably above zero 
at all fields. The general features of the RPM and CPM 
plots are related to the features of the hysteresis loops 
as seen in figure 0] systems with a cliff in the hysteresis 
loop have much smaller RPM and CPM than systems 
with more standard hysteresis loops. 

To more clearly illustrate the relation between disor- 
der and memory, figure ^1 contains a plot of RPM and 
CPM versus "disorder" at the coercive field. "Disorder" 
is quantified by 1/A. There is a definite increase in mem- 
ory as the amount of disorder is increased. Both RPM 
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FIG. 9: RPM and CPM vs. external field as disorder is in- 
creased. The panels here refer to the same systems as in 
figure 121 and figure^ The solid line represents CPM and the 
dashed line represents RPM. All panels have RPM>CPM ex- 
cept in the lowest disorder system where RPM and CPM are 
essentially equal. As disorder increases, both RPM and CPM 
increase. The difference between RPM and CPM is also more 
evident with a larger amount of disorder. 
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FIG. 10: RPM and CPM vs. "disorder" (1/A) at the coercive 
point. The squares represent CPM and the circles represent 
RPM. The increase in "memory" with disorder is evident. 



and CPM are substantial for systems with A > 4.1 which 
corresponds to the point where the snakes no longer ap- 
pear in the domain patterns and the chfF no longer ap- 
pears in the hysteresis loop. Similarly, the speckle ex- 
periment found that this point is where the amount of 
memory is considerable. 



VI. CONCLUSIONS 

The experiments of ref. [iTl IT^ contain numerous ob- 
servations of the domain patterns, hysteresis loops, and 



memory properties of disordered Co/Pt multilayer thin 
films using a variety of techniques. In this paper, we have 
attempted to understand these results by numerically 
simulating the experimental systems. We have provided 
a plausible explanation for all of the results found and 
in particular counter-intuitive results on memory asym- 
metry upon field reversal. We have shown that vector 
dynamics can be crucial in the memory property of mag- 
netic spin systems, without which our model would not 
be able to explain the experimental results. 

Our simulations contain much of the qualitative fea- 
tures of the domain patterns and hysteresis loops. The 
domain patterns for systems with low disorder for both 
experiment and simulation are labyrinthine mazes at the 
coercive field. Furthermore, we observe snake-like growth 
of the domains which are responsible for the cliffs seen 
in the hysteresis loops in accord with experiment. As 
the amount of disorder in the system increases, irregular 
patches replace the serpentine domain patterns and the 
cliff in the hysteresis loop disappears. Because these fea- 
tures in the simulation resemble the features observed in 
the experiments, we believe the parameters used in the 
simulation are close to the actual parameters. It would 
be interesting to further investigate the growth instabil- 
ity in the low disorder regime which should be similar 
in analysis to the ferrofluid case Il3| but without con- 
servation of mass. It would also be interesting to see to 
what extent the analysis of the dendrite problem can be 
carried over to this case 0, . 

We have found that how thermal noise affects the dif- 
ferent domain structures determines how well the system 
"remembers." At finite temperature, the low disorder 
systems become uncorrelated much more easily than the 
high disorder systems; disorder tends to pin domains, 
thereby enabling the system to remember its past. The 
covariance behavior of RPM as disorder is varied in the 
simulations agree with much of the results from the ex- 
periments. 

We are able to explain the seemingly paradoxical ex- 
perimental result that complementary points appear to 
"forget" more than return points despite being governed 
by a Hamiltonian that is invariant upon s —>■ — s and 
Be — > —Bg. The non-invariance of the LLG dynamical 
equation provides a natural explanation for this unex- 
pected behavior. Spin precession reveals itself by decreas- 
ing the correlation between opposite legs. Because of the 
importance of precession for certain physical phenomena, 
scalar theories are inadequate even for highly anisotropic 
materials. From these simulation results, we show that 
the dynamic mechanism is able to explain, at least qual- 
itatively, the observations from the speckle experiments. 
To make this more quantitative, a better experimental 
understanding of the sputtered films is needed. 

Prior to this work, Jagla 0, 0| produced simulations 
with domain patterns and hysteresis loops remarkably 
similar to the experimental ones. In the first of his two 
papers 'l3l , Jagla used a long range scalar cj)"^ model and 
obtained domain patterns very similar to a large number 
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of different experiments. A (j)'^ theory however is not able 
to reproduce major hysteresis loops correctly because (fi 
grows indefinitely with applied field. 

In his second paper, Jagla used a modified model 
that did not allow the indefinite growth of and there- 
fore can produce hysteresis curves that saturate. His 
model Hamiltonian contained all the terms used here: 
dipolar interaction, local elastic energy, perpendicular 
anisotropy, and an external field term. And like our 
work, disorder was introduced in the anisotropy. How- 
ever, there are two major differences between Jagla's 
model and ours. The first difference is that Jagla's model 
in ^3 was at zero temperature and therefore could not 
make predictions on the temperature or disorder depen- 
dence of RPM and CPM. The second more important 
difference is Jagla's use of a scalar model versus our vec- 
tor model. Obviously, a scalar model cannot have preces- 
sion, therefore the mechanism we propose here to explain 
RPM > CPM cannot be applied. And since there are no 
non-bilinear terms in this model (such as a random field 
term), there is no mechanism for RPM to be unequal to 
CPM. However, the scalar model does produce domain 
patterns and hysteresis loops similar to experiments. 

Jagla also showed that when a small random field term 
is included in the Hamiltonian of his original 0^ scalar 
model, instead of a random anisotropy, RPM is greater 



than CPM [13| and the scalar model succeeds in explain- 
ing many of the features of the speckle data even with 
a relatively small field. But an important question re- 
mains: whether or not there exists a scalar theory that 
is capable of explaining simultaneously all of the features 
of the experiment that we are able to do with our vector 
model. 

If our explanation turns out to be correct, this also has 
strong implications for theory, which has often ignored 
the vector nature of the dynamics and used scalar theo- 
ries such as Ising models and cj)'^ theories to understand 
these kind of systems. This opens up the possibility that 
there are other unexplored consequences of this lack of 
symmetry of the dynamical equation. 
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